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I. The Method of E. Zateski. 


fore ere tion, 


The calculation of results of agricultural experiments still repre- 
sents a very important problem, whith is waiting for its solution. All 
solutions, published until now, are not quite satisfactory though often 
they are very useful in practice. The purpose of the present series of 
papers is to consider the two known methods of calculation which we 
believe to be the best and to propose a new one. The order of matters 
in the proposed three papers is the following. In the first present pa- 
per we consider in short the method of ,Student* published in 
»Biometrika* Vol XV p. 271-293 and — a little more fully—that of 
E. Zateski. As the mathematical theory of this method has not yet 
been published, we do so now and show the hypotheses under which 
the application of the same is correct. The second paper will contain 
the theory of the new method which we shall call the method of para- 
bolic curves. In the last paper we shall publish some tables which 
are necessary, when the method of parabolic curves is applied and also 
several examples with applications of the three methods, which will 
illustrate their different advantages. 


I. General Considerations. 


The first question to be considered here is: how to verify a me- 
thod of calculation? Or, how this method ought to be in order that we 
may find it satisfactory. The answer is: the satisfactory method of 
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calculation must 1) eliminate as far as possib!e the inevitable errors of 
experiments and 2) must be simple enough for frequent use. Certainly 
there are no methods which will eliminate all errors of the expe- 
riments. But knowing some of their sources and the characters of the 
resulting errors, we can adapt a method which will eliminate just them. 
It is Known that the main sources of errors in agricultural experimen- 
tation are: 1) inequality of the soil used for experiments. and 2) the 
inevitable random error, due to a large number of sources, each of them 
giving only a very small error. Only the first of these two sources of errors, 
is interesting, because, if the action of this source is eliminated, 
increasing the number of experiments we shall easily eliminate the 
other. Thus we shall ask whether and how a method can help us in 
eliminating the errors of experiments due to the inequality of the soil. 


We must confess that the expression we are using, namely the 
,inequality of the soil* is rather vague, and we do not precisely know 
when the soi! is to be considered not equal. We know only that very 
often the yields of different sorts of cereals compared in a field expe- 
riment seem to increase more or less systematically when we approach 
a certain point on the field. In these cases we are inclined to think. 
that the quality of the soil on the field is changing and that the men- 
tioned increase of yields is-due not to random causes, but to changes 
in the soil. This assumption is confirmed by the so called blind expe- 
riments, when all plots of the field are sown by a single sort of cereals. 


The question is now to find out what the changes of the soil are. 
In other words: let us assume that a certain sort of cereals A be sown 
“1 plots 


PRO pet aia ee Se » An (1) 


situated somewhere on the field, and that the respective yields were 


Diy Yop ee ee, Yn; (2) 


assuming further that we know the yields of remaining varieties on 
other plots, we ask the question: Which would be the yields of the 
variety A if it were sown on all other plots of the field besides (1)? 


We have no possibility to answer this question accurately and we 
only assume that if f. i. round a certain point P» the other compared 
varieties, say B, C,...., Dare giving yields which generally are larger, 
than the yields of the’same varieties on more distant plots, then pro- 
probably the yields of the variety A on plots round P would be also 
Jarger than the yields on other plots. But what would be the mentio- 
ned increase of yields of A? Would it be proportional to the increase. 
of yields of other varieties? We do not know: the experimental data 
are too confused owing to the random error. In some cases they seem 
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to contradict all hypotheses of that sort. Nevertheless all methods of 
-calculation of experiments assume that there exists a certain level of 
fertility of the soil, common to all varieties and, if we could know the 
true yields (not affected by random errors) of two different varieties on 
-all plots of the field, the difference between respective yields of ihe 
same plot would be constant for all plots of the field. 


AS we are now selecting and comparing the varieties which are 
very-much alike, this hypothesis is often very nearly true. The practice 
shows that we may accept it almost always, when comparing the cha- 
‘racter of varieties measured by some ratio, for instance when the cha- 
‘racter under consideration is the percentage of sugar in beets and others. 


But the hypothesis that the difference between the true yields of 
different varieties on all plots of the field is constant or approximately 
$o, is not sufficient for the calculation of changes of the level of soil- 
fertility: we must accept some further hypothesis concerning the man- 
mer of changes in the soil. There are many different kinds of such 
hypothesis and each has a corresponding special method of calculation, 
-so that accepting or rejecting the hypothesis we must do the same 
with the corresponding method. We shall start by considering briefly 
»student’s* method and the hypothesis on which it is based. 


II. ,Student’s“ method. 


»students* method assumes that the field may be divided 
‘into parts of the same size, on which the level of fertility is sensibily 
‘constant. Assume for simplicity that we compare only four varieties 
A, B, C, D and that they are sown on 20 oblong rectangular plots in 
one row in the follwing order 


Serer eee? C.D. . A, 3B, 6, D 


-so that each variety is repeated on five plots. ,Student* proposes 
in such cases to divide the field into five equal parts, which we shall call 
‘series, so that on every series every variety would be sown on one 
single plot. He represents further every true yield of a variety by 
-a sum of two components: the first one X; is the level of fertility, 
‘common to all varieties but varying when we pass from one series to 
another, and the second Xz, characteristic for the variety, constant for 
all yields of the same, independently from the series on which the 
‘variety is sown. So for instance the five true yields of the variety A 
‘would be 


X, +X A, Xo KA, Xe tXa, Xe PX A, XX A, (3) 
ithe five yields of the variety B: 
X,+Xep, X»-XB, Sg7A B, XsXB, X;+XB, (4) 
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and so on. The difference between respective observed yields of every 
vuriety and the hypothetical true yields of the form (3) and (4) ,St us 
dent‘ ascribes to the random error. 

The problem of calculation consists in estimating the four con- 
stants Xa, XB, Xc, Xp, together with their probable or mean errors. 
The solution of this problem is perfect and the statistics introduced by 
»otudent* for estimation of X’,s and their mean errors are easy 
to calculate in practice. 

In order to judge whether ,Students* method is quite satis— 
factory we must ask how far it eliminates the differences in the soil. 
As mentioned we do not precisely know, how these differences are,. 
but we are rather inelined to think, that the fertility of the soil is 
varying in a, say, more continuous way, than it is assumed by ,,Stu- 
dent‘. The division of the field into several series of plots on which 
the fertility is assumed to be constant, seems to be artificial. In our 
example the fertility of thefourth plot is assumed to be the same as 
that of the first one, and different from the fertility of the fifth. We 
would rather suppose, and this hypothesis can be easily confirmed by 
observation, that the difference between the fertility of the first plot. 
and that of the fourth is generally larger than the respective difference 
of the fourth and the fifth. This objection becomes obvious when we 
try to represent the fertility of different plots by a line in such a way 
that the abscissa will be the number of the plot and the ordinate — the 
respective fertility. If the ,Students* hypothesis were true 
then the fertility of a field would change as shown on the following: 
drawing I. 


Considering this drawing (which corresponds to an actual experi- 
ment) we are pretty sure that the fertility of the soil is changing con-- 
tinuously, rising at the beginning and decreasing at the end of the 
field. We see also that the application of ,Student s* method 
would lead to an error in the estimation of single yields of varieties. 
A and D. 

Of course, ,Student* knew the mentioned difficulty, but. the 
scheme of experiments he dealt with did not allow him to avoid it with uti 
introducing terrible complications in the resulting formulae. We must. 
not forget that all assumption concerning the character of changes in. 
soil will give us only an approximation to the facts of observation. 

The ,Student’s* method gives just only an (often very good) 
approximation to the actual changes in soil-fertility, but as these seem. 
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to be rather continuous, it is quite natural that the attempts are made 
to find a continuous curve which, fitted by a certain method to the 
experimental data, would represent the theoretical level of fertility of 
the soil. Such methods have been invented almost at the same time 
by several scientists from which we shall mention FE. Zateski and 
mins cheriich. 

The methods are very simiJar so we shall consider only one 
‘of them namely that of E. Zateski. We shall see that from certain 
points of view it is also not quite satistactory, though it often gives 
‘very good results when applied in practice. 


Il. E.. Zateski’s Method. 


a) General Description. Neither ofthe two authors occu- 
pies himself with mathematical hypotheses concerning the character of 
‘changes in the soil, which would lead him to the proposed methods. 
They simply follow the good common sense and state that if we calcu- 
late the mean value of several yields from neighbourixg plots we shall 
partly eliminate the random errors and get a number, which will be 
often very near to the mean of yields of the same varieties which 
we would get cultivating them in the same condi- 
fiemes on whe plotiin the middle of the group under 
consideration. Thus the described mean of several, say neigh- 
bouring yields would give us the required level. To complete the 
description of the method of E. Zateski we notice only that to have 
a level common for compared varieties we must use them all for cal- 
- culation of the level in every point of the field. 

The position is quite clear if the number of the varielies is odd. 
Assume for instance that we compare five varieties A, B, C, D, HE sown 
as usual in the order 


A, B,, C, Dy, &,, A:, B,, C., D,, Ag. (5) 


‘The subscripts indicate here the order of. the plot sown by the same 
variety. We have no means to calculate the level for the first two plots. 
‘The level for the following will be given by means of formulae 


ee he Pe DBAs, Ce Der Birr ayy Be uh) 
5 5 P 


ff the number of varieties is even, for instance. if we compare four 
varieties A, B, C, D, we calculate the level by means of the formulae: 


Ye AE Bi tC, +D,+'/2 A, "ls BerO TD rArr'/:Bo, ete. (7) 
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We. see that in both cases we have no levels for the first two and: 
the two last plots of the field. Assume now that the level, say Z has 
been calculated for all plots sown by the same variety A. Taking the 
differences between the observed yields of A and the corresponding | 
levels ” 


Vai =Ai- —' Li (i Oras) A + | lll 


we have the postitive or negative excess Vai of the yield of A upon the 
level. Clearly Vai is very much affected by the random error in A, so 
we try to eliminate this error in calculating the mean value of all Vaiss: 


nw 
18 Cee YO i (9) 
%. 


where * — 1 is the number of plots sown by the variety A, for which 
it is possible to calculate the level. If we wish now to compare the: 
productiveness of two varieties say A and B, we take the difference 
Va — Vs. , 
b) The Fundamental Hypothesis. For simplicity we 

shall assume that the number of varieties we compare is odd 2 v + 1, 
Assume further that the varieties are identical and that we are dealing 
with true yields, or that all random errors are eliminated. Then applying 
the method of E. Zateski and calculating the level we should find for 
every plot the number, which measures its observed (and in conformity. 
with our hypothesis, true) yield. We shall see that this does not 
happen in all cases of practice. Let us call the observed true yields of 
our (identical) varieties from subsequent plots 


28 Cee @ a Me C (10) 
The level for the i—th plot would be given by 


j Y ; 
Li = MS : 
v+1 ey ; Yitk (11): | 


Stating that it is equal to the observed true yield Yi 


1 i 
Yee v (19%: 
Seco re a (12) | 


we get the condition which must satisfy the fertility of 
the soil in those cases when the application of Za- 
Jeskis method gives no systematic.error, 

The condition (12) is a homogeneous equation in finite differences 
with constant coefficients and can easily be solved. For instance if v be 
equal to 2, the general solution is given by 
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Ye: =a-- O-- eri? dry* | (13) 
where 
le AT 
Yeo og CCN (14) 
2 
: as 
ee wens Coe (15) 
2? 4 


and the coefficients a, b, c, d are arbitrary. Generally the solution of 
(12) is given by 


Ve Sav beer* Pore Py. . (16) 
where a and 0 are arbitrary, the r’s are the different roots of the equation 
1+2z+2?+....+ 2¥-1 — 2, zv + 2vtl gvt2 +... zv = 0 (17) 


and the P’s — the arbitrary polynomials in x of the order less by unity 
than the multiplicity of the respective root r. 
We see that putting in (16) instead of P’s zeros, we get the equation 


Yx = ax-+b (18) 


of the straight line—thus, if the true level of fertility is represented by 
a straight line the method of E. Zateski will not introduce any syste- 
matic error. The same will happen in much more general cases, namely 

when the P’s are not equal to nought. To have an idea about the shape 
of such curve we have drawn it for the case v=2. The arbitrary coeffi- 
cients a, b, c, d are chosen in order to allow the curve to pass through 
four given points, represented on the drawing by larger dots. 


Pig L /\\ 


Fig. 2 


We see that in the middle of the drawing the curve is pretty 
smooth but at the ends it is rather rough. The bounds at the ends of 
the curve are due to the circumstance that r, and r, are negative and 
their powers r* are changing the sign when we pass from an even x to 
the following odd one. 

x 
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What we have seen in this special case may happen in the general. 
If the curve corresponding to the equation in finite differences (12) is 
not a straight line, then the deviations from linearity or bounds may 
indefinitely increase in both directions. In fact let us put in (12) i=xtv 
and write the same equation for three consecutive points. We shall have 


2Qy 


Yuty Sige x ry (19) 
Qv+1 

Yeo em : a Yetk (20) 
1 2yte 

Ye4y42 = pune ‘ 2 : Yutk <2}) 


Subtracting the second equation from the third one and the first from 
the second we have 


1 
Aviva = VYeovie — Yet = Cia (Yx+2y+2 — Yx41) (22) 
/ 7 1 
an = YetyH— Yety= Ova. (Yx+a+1— Yu) (23) 


Subtracting again the last equation from the preceeding one we have 


AQ) = LO (Yesoned — Vedi = (Yee (24) 


By tt 2y+1 
Or 


Yu+toyt2 — Yatoyti= Yeti— Fx + (2441) ARM x+y (25) 


If the so-called first differences Aa=Ya+i1—Ya are for every a constant, 
the points with coordinates (a, Ya) are on the straight line. Suppose 
now that the ordinates Yx+y, Yxtvti, Yxty+2 do not  corres-~ 
pond to points lying on a straight line, so that two first diffe- 
rences in brackest on the right hand side in (22) (23) are not equal. Let 
their difference be equal A‘), —q>o Then we shall have 

Yu4-2+2 — Yutoy—i = Yoti— Ya + q (2H) (26) 
or with other words—the difference }x+2 v+2 — Yx+2y+41 will be larger 
than Yxti—Yx by a quantity 2,11 times greater than g. The position 


s clearly illustrated on the drawing. 
Asriver 
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Here the points Ax, Ax+ti, Ax+y, Axt+y+i1, Ax+y4+2 are drawn to be not 
very far from a straight line, but while the differences of the ordinates 
of the two first and the two last points of this series are both equal to 
‘/,, the difference between ordinates of Ax+y and Ax+ y41 is supposed to 
be zero. Now, as it follows from (26) whatever the ordinate Yx+te2y+1i 
of Ax+2,+1 is, the ordinate of the following point must be larger by y+1 
or, if, as in the considered case y=5, by 6. 


In this way we see that if in a certain region the first differences 
Ax+, for consecutive values of X are varying. even very slightly, the 
variation of ordinates in points Ax+2,y will be considerably altered as 
compared with Ax, Mutatis mutandis the same can be proved for the 
case when the number of varieties is even. 

This property of the considered curves is not very convenient, as 
we are rather prepared to believe that the changes in the soil are smooth 

The second inconvenient property of the curves under considera- 
tion consists in the fact that they do not have points with ordinates 
exceeding ordinates of y neighbouring points on every side of them. 
In fact if it were | 


Yi-k < Yi; Viner < Yi(k = 1,2...) (27) - 


we should have 
1 


ye Yer (28) 
k=—y-+1 


mM 


4 ¥ 1 [y 
Leg 2v+1 Srey ben Ae 4y Vike Vita 


and the calculation of the level by E. Zateski’s ,method would 
introduce an error. In the same way no point of the curve can have an 
ordinate which would be less than the ordinates of the same neighbouring 
points. If we turn now to practical data, we shall state that very often 
we meet smooth levels of fertility with two or three maxima and minima 
on the field containing some 5(2v+1) or more plots. We see that in these 
conditions no curve of the considered family can represent rigorously 
the level of fertility, as these curves have either no maxima and minima, 
or much more often than it is on the field. We find now the source of 
systematic errors in applying the method of E. Zateski or that of 
Mitscherlich: If the true level of fertility has in a certain point 
a maximum, such maximum is by using the considered methods underap- 
preciated and thus the yield of the variety sown on the same plot su- 
perappreciatied. 


c) The Estimate ofthe Standard Error. 


Calculating the estimate of standard errors of results of calulations 
by the method of E. Zateski we shall use several general formulae. Let 
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Wl, Us, ere a un ; (29) 
be some variates with equal standard deviations g and 
u=)\ Uy + he Uo t....+dnun (30) 


a linear function of them with constant coefficients }.. It is known that 
if the variates (29) are not correlated, the standard deviation of their 
function (80) is given by the formula: 


+ | 
Gu = / x! hi (31) 
i=l 


The second theorem is the following. Let us consider different tncor- 
related variates 


i gl gy Rene STL ae (32) 
having in a population different means, respectively 
Ui, 4p) |e OE ne ae (33) 


and the same standard deviation co. Assume further that the set (82) of 
variates can be divided in groups each containing only different variates: 


a a HENS | _ (34) 


but where not necessarily all members of one group differ from these of 
another group, and such that for different groups a certain linear function 


U= Mt ‘Ui, -- hs Ui; + eeee + An Win (35) 


with known constant coefficients jn has the same constant value. Now 
we shall consider a variate defined by the formula. 


U aaa Di pahe Uy ee oa a eae Oe (36) 


in which, in order to get the values of y it is necessary to substitute 
the values of any group (34) instead of v’s. Clearly the variate y.has 


its mean equal to y, and owing to the hypothesis that all w’s have the 
same standard deviation 5 — a standard deviation 


: piey nt aah | 
5) = // Ye (37) 
b=1 
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Suppose now we have a system of observed values of m groups (34) of 
u’s. Putting them into (36) we obtain a system of values of y: 


VU iL U2, Core ia) Um (38) 
where 
n 
k=1 
If Uk means 
m1 
U, =— 2 Mir (40) 


m i= 


then the mean of the sample (88) will be 


Ret 


Now consider the deviation of U; from its mean in the sample. We have 
Cf Voss Une ee Ciao) 


or, owing to the property of (35) 


: n = n ee - 
Ce EW Beds Ap Ree) | Vee pecan U (42) 
k= k= 
Now we substitute instead of Uz its value (40) and instead of y its 
value, which may be written 
= 1 nt ze as 
Uy Sa Dy, S Ae YR (43) 
P= 10k =f : 


Tl 
In fact, owing to our hypothesis ) )e ui has the value y independen- 
k=1 
tly from the subscript i. Thus we have 
rca pet ml ia] a 
: i telpscs 
MU—V= YX ft ee ng| OPES: Ge oe te | (44) 


k=1 j=l. k=1 


The differences like U,— V serve us for calculation of the standard de- 
viation in the observed sample of values of y. Let us calculate the 
mean of (Ui — V)? in repeated samples (38). We notice that the diffe- 
rence vy, — V is a linear function of all observed values wir. Owing to 


our hypothesis, not all variates in one group (34) are necessarily 
7 $s 
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different from the variates in the other. Therefore without further 
assumptions we cannot write explicitely the coefficients of y— V 
expressed as a linear function of different w’s. Nevertheless we may write 


Uy —- vy = 32% (a — v) (45) 


whare @s are certain constant cefficients, which can be determined 
when ten character of groups (34) is specified and which may vary with 
the subscript i. The sum in (45) is extended over all observed diffe- 
rent values of uw’s. Now the mean value of (yi— V)? can be easily cal- 
culated, the variates ui being independent Denoting the mean value 


of every variate X by {x} we have 


Ke — vy} = ek — WP 4 22aa(u—w (ua) } 


‘where the second sum is extended over all products of differient @ (u—u). 
As 


(a a) ty = 3 (47) 
{((a — u) (u— u)}=o (48) 

the u’s being uncorrelated we have 
{( — yy} = date (9) 


where the sum Y¢? can be dependent upon the value of 7. We see that 


1 
ee (vi — V)* may be considered as an estimate of co’. To imovrove the 
tod I, 


accuracy we may calculate the mean value of all such expressions in the 
sample (38), say 


limp srr nes (50) 


It may happen however (it does for instance when there are no common 
terms in groups (34)), that La? is independent from the subscript i. Then 
the expression (50) is nothing but the ratio of the ordinary variance of 
v’s in the sample over Yv2. In any case we shall retain the above 
notation and write 


Estimate of 52 = a’? (51) 


This estimate may be used for the calculation of the variance in 
reneated samples or in the population of y and other variates connected 
with vy. For instance the estimate of the variance of y in the population 
~ will be 

12 
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Estimate of os — 2 2 6 (52) 
Uv 


The formulae (51) and (52) will be further of frequent use. The meaning 
of (51) is as follows. 1f a variate yis connected with the 
moeurrelated variates ua’s as ‘specified, and if u's 
have a common standard deviation oa, we get the esti- 


mate of g2 in calculating the mean value of 


eD 
the «’s being the coefficients of different uw’sinthe 
expression of yi — V. An example will clear the matter. Suppcse 
all variates uik be identical with a single variate w and that n=1: every 
group contains only one element: the value of u observed in one trial. 
Suppose further that }, = 1, then y =u and yi will be the value of w 
observed in the i—th trial. 

What will be the a’s? 

Turning to (44), putting n=1 and omitting the subscript k everyw- 
here, we have 


ED 


m1 (ui i Pe un) Cae 


nt 


X'(uj — a) (53) 


Vik Ness 


where the sum )’ is extended over all values j/ = 1, 2.... m except j =i. 


' m—1 Or i. 
Hence one of ’s is equal to and the remaining ones to >, indepen- 


~ 


dently from the value of i and so 
ee ee (54) 


Thus the estimate of the standard deviation of the variate w is 


le = 
oi ig a : e 


the very well’ known result. 


Now let us turn tothe method of calculating of E. Zateski 
Suppose first that we compare 2y+1 varieties, the number‘of plots sown 
by each of 2y varieties being y and the same number for the last ,stan- 
dard* variety 4+1. The yields of the ¢-—th variety we shall deno!e by 


Hl tgy ht: Bs towne be : (56) 


the yields of the standard variety — by 


VOUS YQ 2 ned X21 (57) 
| 1s 
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—— 


where the second subscript means the order of subsequent plots sown by 
respective variety, numerated from one end of the field to the other. 


Ift <y, x > 1 then the level for the x — th plot sown by 7— th 
variety is, as we have seen, estimated by the mean 


1 2y t+ by \ x! 
ASRS ) Yex—1 + + DVR x J (58) 
Re: k=t+y+1 k=0 


The hypothesis accepted by this method is that the true yields of every 
t{—th variety exceed the level on respective plots by a constant quantity 
which we shall denote by Uz, with the subscript ¢—correspondent to the 
variety under consideration. Thus the difference between the observed 
yield of the ‘—th variety and the estimate at the level 


Qy t+y . 
Oe pal Gh een k py Yee ee Yr.« (59) 
yO N pty k=0 


is considered as a particular value of a variate Ut with the mean Vj, x 
The expression of Us, may be written also in the following way: 


he a ee 
ie ot x Yr, fh cae re % Yr«x + Yr,x (60) 
a k=t+y+1 k=0 k=t+41 
or simply: 
2y al > 
Of eo = —— — “<“(f, x j 
vs AE ee Spe ee | (61) 


“V+, x can be calculated not for every variety and every plot. Ilft< y then 
such an exceptional plot, for which the level and Ut,~ cannot be calcu- 
lated, is the first plot sown by the ¢—th variety. If /~v--+1 the same 
is for the last plot. For the standard variety the variate Ut,x can be cal- 
culated for all plots except two — the first and the last one. The varie- 
ties the y— th and the y + 1 — st one have the Ui, for all plots sown 
by them. All this can be easily seen on the drawing. Assume first that 
o< t< y and find the expression which is calculated for the estimation 


of Uz, It is the mean value of all Uz, x(x = 2, 3.... x), which we shall 
denote by Vi. We have 
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3 ¥. 
y= NM Ytx — 
(Sue (hl) Soe 6 
62 
; U 2, 4 gh, (62) 
meas iy) = Garg es eae ae ee ge 
y ‘ x2 Shetty k=0 k=t-+1 
Further in the analoguous way: 
g Me 
dy x \ % yl 2y 
yx i Ni ibaa, eo cea AS sy 4x Tig c rE 
bee Qyliy. Mon. (Qyth, 7 u Yat X Va, (63) 
tet MPA vovapioe | SNS Baas () k=yt1 
2Qy % ; 1 % ¥ Yan 
Vy — ee YH, 2» eg ete oe » by Yrix + > Yr, x + Yo, a+ 
(2v+1)z {2y-+1)y, 
and at last, for v4+1<t (64) 
. oy %-1 
4 = —-————_ 8 Yy. oo 
(2y+1) (¥—!) Pe ue 
35 
' v1 / t-1 ai eS (65) 
Vv va a aT r 
Be iy ee ee Ya 
i CP iecel eat k=t+1 k=0 


In order to use our previous results, namely the formulae (51) and 
(52) we notice that E. Zateski’s method of calculation consists 
really in dividing all variates Yz,x into groups, for which a certain 
linear function with constant coefficiénts, namely the function (60) Vz,x« 
remains constant provided we put in it instead of Yi. their mean 
values. In fact we remember that one of our hypotheses was that the 
excess of the yield of every variety over the level is on every plot 
constant and characteristic for this variety. 

Consider for instance the »-th variety. The corresponding groups 
(34) are the yields used for calculation of levels for plots sown by the 
same variety, that is to say 


Yan yi Vine odes Ne hee 
Ypres Yisge Yoroe vee Yo yes 


(66) 


Vera. tee Long ecns Loy 


We notice, and that is true not only for the v-th variety, but in 
general, that different groups (66) have no common terms. If it were 
the ,true* yields and if the fundamental hypothesis of E. Zateski’s 
method were fulfilled, the mean values of these varieties would keep 
constant the function 
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ee! 2y a i 1 y—1 me 2y a. 
Uy x = oy +1 Yy,x tg oy +1 | ae Yr,x + ier Yr,x (67) 
= =y 


analogaus to (35) independently from the value of x, where Yi,x means 
the mean value of Yi,x. Hence, as we suppose that there is no 
dependency among errors committed in measuring Y’s, that is to say that 
in repeated samples there is no correlation between Yi,x — Y+,x and 
Yz,& — Yz,£, to have an estimate of the common variance of Y’s we 
have only to calculate the mean value of all terms of the form 


Cur, —Vt y? 


ya? 


(68) 


the a’s being the coefficients of different Y’s, in the expression of 
Ut x — Vt. As the groups (66) dont have any common terms, Xa? 
has a constant value independent from the value of x, but changing 
with the subscript 7. 

Let us calculate \g2 for the two cases, when ¢<y or y4+1<? 
and when t = y or ¢=y44 1. Suppose first ¢ < y. 


We have 
2 (4-2) 2y aed - 
Vixe—Vix Yt, «— ——— Y, Yin) — 
oe LOA) Gy ME (Eee (2, ak SF out e 
4-2 | 2y Luho ob i+y 
pal Ye, ca See Yas) + 
F tae dane) ’ a ’ ’ 
CD Ca-t) \ pay k=0 kot (68) 


{ x—1 % \ ey t—1 f+y 
oe ("s+ ¥ }( Y Yire41+ 2.Yi2+. 3 Yin), 


(2741) (4-1) \ 5 koxti/ \ist+y+1 i=0 i=t+i 


ie me x 

where & -+ % indicates the summation over all values of k=2, 3... % 
Rees 2 k=x+1 

except for k =x. Then we have in this expression no terms with equal — 


Y’s, thus the coefficients are the @’s required. There are theiefore 


. Ao Eh 2y (* — 2) 
1 term with coefficient f, equal to (v4) Gad 
29 
#——2 terms -, ” “ go, XO) ogee 
Par 
2yv » ” » 9 9 ~ Oy+i@—D 
. ib 
2y (x—2) ” ” ” ” ” » (2y+1) (4-4) 
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Hence 
4y? (x — 2)* (x— 2)? 4y2 
Vy2 = aN Fy ee 
Ee (2y+1)2 Git : (2y+4)2 (4—1)2 peer? (Qy--1)S (42)? 
(69) 
2y (x%—2) 
DVttae se = = 
sig ) (pyar ieee e-1) 5 (y+) C—) 
and the estimate for the variance s* is thus 
4 2y+1) (y—1 
Estimate so? = iba lads Ms (vVix— Vt)? (70) 


2y (%—2) 


Evidently there is no essential difference in reasonings or in calcu- 
lations when ¢ is supposed to be = Y. The only difference is that if 


t = VY or ¢ = V+ 1 the number of plots for which vz,x can be calculated 
is not x — 1 but x. Thus o* may be estimated by means of the last for- 
mula for every value of ¢ except for = ¥,¢ = V+-1. 
When applying these two varietes the ¥Y — th and the v-+ 1 st we 
may use the following formula 
(24y+) x 


E ti te se oS A) pee 
stima oe (Vt,x Vt) (71 


So we obtain (2 ¥ — 1) («x —1)-+ 2 ~ different estimates for os. We can 
improve the accuracy of estimates considering as the estimate of 5° the 
mean cf all just deducted estimates. 

1 . (2y-+41) @%—1) 
(2y—-1) (x—1)+2x 2y (x—2) 


(vt«%— Ve)?+ 


Estimate s? = 


(2y+1) » : x 
a (Us,x —Vs)* 


Vv 
oad 
ay C4— 


Where the first sum is extended over v¢,x corresponding to all varieties 
except the V — th and the V + 1—st and the last—over the remaining 
ones. 

The found expression may be-used for estimation of the standard 
error of V’s and others. 

When comparing varieties we are mainly interested in differences 
Vitp — Vt which show how far the i- p-th variety is better or worse 
than another the ¢-th. Therefore we sha!l use the previous results for 
calculation of the standard errors of such differences. The method will 
be as indicated above: We shall express Vitp—Vzi as a linear function 
of Y’s, calculate the coefficients 4 of different Y’s in it, square them 
and multiply the sum of squares by the estimate of o*-the Tesla will 
be the estimate of the variance of Vttp — Vit. 
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We must distinguish several special cases. Some of them are 
arising clearly from the differences of expressions of Vz for different 
values of ¢. Besides two special forms of the variance for the same 
expressions of V’s will arise due to the special mutual position of com- 
pared verieties. Let us consider the difference Vt+p — Vi. 


Case J]. We assume: ogc i~t+p<y. 


We have 
2y % 1 % 2y 
Vitp — Vi = SONY £25 oe ae ee bs Yixr—i 
(vt) G—-D ._o (Cyt) GD) 4 Ge 
i-+p-1 t+pty Qy % 
+) OY Yt Yia|— ee (73) 
i=0 fated (Vt) @—-D 
1 % 2y i—1 ity 
+ $y > Yi, x—1 + > Yi,x t >; Yi,x 
Coy (Oar ale Vee ame i=0 $=t1 1 


Considering this expression we see that a) its second term has 
members with Yz,x.’s, b) its fourth term contains members with Yt#+p,«’s 
and at last c) that the two terms the second and the fourth have seve- 
ral equal members with different signs which will vanish. Our purpose 
is to calculate }’s, that is to say coefficients of Vitp— Vi represented 
as a linear function of different Y’s. Therefore we proceed to put 
together terms having the same Y’s.. As ¢-++>p<v it must be ¢+p<ié+-yv 
and so the fourth term of (73) contains Yt+p,. with the same subscripts 
x= 2,3..% as the first term. Putting them together we have the sum 


2y % 


. A a 
EY, Hae DE he Pr bt 
(2y+1) (x—1) a te P (2y+1) (x—1) ci 


| % 
Yitpis = a di Yt+p,x (74) 


x=? 


In the same way we state that every term of the third member of (73) _ 
has a corresponding term with the same Yz,x in the second. Putting 
them together we cbtain the sum 


Saba (95) 


Consider now the remaining terms of the second and the last members 
of (73). It is-easy to see that for t+ p+v+1 <i < 2v and for 
o<is<ct-+y the subscript x runs in both members over the same 
values. Therefore the corresponding terms, being equal and with dif- 
ferent signs, will vanish. For the remaining values of i, i. e. for 
i+v+1<i<t+»-+p p the values of x in both members are different: 
in the second member x runs over values 2, 3,....% and in the fourth 
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over 1, 2, 3,...%x—1. Consequently in the second member all terms 
‘will vanish except with x =x and in the fourth—all except with x= 2. 
Thus the two considered members of (73) will be reduced to the sum 


1 t+y-+p 


ore ee en Yix— Yi 76 
(ay bl) (ead) hihi t,%—— Yi,1) (76) 


‘In this way ihe difference (73) may be written: 


1 % 1 % 
Vip —Vi=-—- » Vip xo ee le py Yt, x — 
‘ot cee i eS 
(77) 
1 t+y-+p 
ay CANAL) (petite | ro meee 


t=t+y+1 


where there are no terms with equal Y’s. Thus the coefficients of this 
expression are the required 4’s. There are 


%—1 of them equal to 


%—1 
1 
== yt 7) ) ” a ara 
(78) 
i 1 
p ” » ” Saou C1) (een) 
1 
eee? yt 1) Gt) 
“Therefore 
2 2p 
Ea he essere mean (79) 


u—1 (2y+-1)? (4—2)2 


‘Multiplying this expression into the estimate of o? we get the estimate 
of the variance of Vitp — Vi. Clearly the same formula (79) is valid 
ieee tt p << 2y. 

We leave to the reader the care to prove with quite analogic 
ecalculations that in 

Case II a) when ¢~v, v+1<t+p, <p 


Ne et! % 
Vitp—Vi= aay, > VYt+p,<— ma % Yixn—- 
as | K " 4=2 (s0) 
% 1 tty 


we eee, Se Ss > Y; ic Yt 
C2) ey ee ‘ ; 


—_— 
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and thus 
2 2y+1—p 


ye 2 Penta Wiese Dalek: 
ON ARON (2y4+1)2 (y—1)2 


Rech 
The formula (81) seems to be quite different from (79), but we 
shall see that this difference is rather apparent than essential. In fact 
in the Case I p means the «nearest distance between plots sown 
with varieties the ¢-th and the ¢t-++p-th. In Case Il a) p is suppo- 
sed to be >v and as we compare 2v-+1 varieties, the nearest distance 
between the plots sown with the two varieties is not p but 2v-+ 1—p. 
lf we denote this distance by d, the two formulae (79) and (81) may be 
written in the common form 
one 2 2q r 
S42 == —  ——_____ 2), 
det Ti) Gt? Go pe 
which is valid for oxd<v and when the nearest 
plots sown by two compared varieties are not 
separated by the plot sown with the y»y-th variety. 
If the nearest plots with two compared varieties are separated 
by the plot with the v-th variety we have the 
Case IIb. The reader will easily find out that in this case 


‘ Qy 1) 4 1 
—Vix= Y; Src Y; pe aI — 
1 1 nye 2y 
2 Yt OD OV a ee 
Cr) Gay ats eee 
1 | t-1 (83) 
Se ie Yii— Yi, x»)+ 
= l, ’ 
Cd) GD [po aa, 
t--p—1 tty . 
+ YY Yn-—-Yiy+ %L (C¥21—Yi,x) 
i=tH i=it+pH 
and therefore . 
: 2 2y+1+-p 
pa — (84) 


%—1 (2y-+1)2 (y—1)? 


This formula is valid when p<v and when the two compared 
varieties are separated by: ithe v-th variety. We see that the compari- 
son between these varieties is more accurate than in preceding cases. 
This is due to the fact that the difference (83) is depending upon all x 
yields of the two. varieties while in the Case IT and Case Ila 
only »—1 yields have been used. 


), ) 
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Deducing the formulae (82) and (84) we have assumed that not 
one of the compared varieties is the v-th or the v+1 stone. 


owe Oil) fev, p = 1, 
We have easily: 
1 % J 


Yyiy Pe om ” 2 Yy,x are Qy+Dx (Yo, x4] =~¥'o,1) (85) 


s==| ine 


| ] 

B= 2{—+ 

4 | * . y+ 1)? x? | fF) 

Case IV. The reader will easily persuade himself that when 

the y-th or the v+1-st variety is compared with any other /¢-th variety 
we have: 


] y 
PSI a ee ene pecs pat 8 
f cre Qy+)x aay CoH) 
if << v, we have for the difference Vi — V, the following expression: 
y. y | 2yxtx—| es 2y 
eee rt Yh ee re ee Ye pe — 
Wet fey GoD 2. Ot oe 
Ryd — T)-F is 1 2y t—I 
ae . nora Yy ox me p> MAP wt y Yin a (88) 
@yTDe G—) 5 Qy+!)x yy 1S ; 
y-1 i+y Ve] % ay 
ot Yi. tive | BV, | Senea enn py y Yio2-1-+ 
i=t+-1 t=y+1 aye Onl) 2.5 i=t+yt+1 
t—1 y—l ity 
+ py Yi,x+ > Yi,x+ x Yr 
i=0 i=t+1 i=y+1 


In other circumstances contained in Case IV this expression 
is quite analogous. 

We shall now occupy ourselves with the case when the number 
of compared varieties is even 2v. As we shall see this case is more 
complicated than the previous. We start with noticing that the groups 
(34) are now overlapping as for instance for the v-th variety they are 
as follows: ; 

FON aY 15.0 Loy tLe 0.2 
VO 1 ey Sed eh 08 


Fe tees) ee rn (89) 


YO, Yijny +0» Y2y—1,x, YO, x+1 
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where again the first subscript 1=0, 1, 2... 2v—1 corresponds to the 
variety and the second indicates the order of the plot sown by the same. 


We shall start with calculating the estimate of the variance 3? 
common for all Y’s. Supposing ¢~v we shall have the expression for 


the level 


1 1 : 2y—1 i+y—1 | : 
Lies zi b Yit+y,y—1t >» Yixn-1+ YS Vix oy Ytty, | (90): 
E i=ttyt1 i=0 


and for the excess of the yield of ¢-th variety on the x-th plot: 


et 1 2y—1 i—| 
Di aT Yea "(4 Yee © Yix—-1+ % Ynaxt 


t+y—l 
» Yi, x + 4 hae 1 
i=t-+1 


As vx can be calculated for x=2, 3... %, to have the meani 
excess Vi we have to sum (91) over all these values of x and to divide 
the result over » —1. We notice that owing to the fact that the groups 
(89) are now overlapping every two expressions vzx« and VU¢.«ii1 have 


1 
a common term, namely — [ay Yi) es 


Therefore, we have 


y Q2v—1 bs ines 1 1 y y. 
t 2) (wt) ve t,x Dy (wae (Ytty,1+ ttyy) 
(92) 
x—1 2yv—1 % t+y—1 % 
+ YY Yoyext YF >. Yge— te ae Yas 
t= i=t+y+1 x=2 i=0 22 


b 
where the sum ¥ (¢) is extended over all values i= a,a+1.... b, except i= t. 
i—@ 
To have the estimate of the variance of Y’s we have to follow 
the formula (51) and to calculate the mean value of 


COL haa ed. 


Ya? ee) 


the «’s being the coefficients of different Y’s in the expression of Uix—Vz, 
It is easy to see that Meis now dependent on the value of x. We 


have in fact: 
ae 
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(2—1) (x—2) | feat le a a 
U12— Vi= 2 er eae | x Yt } WMYi2|\— 
2Qy (x—1) 2y (x—1) t=t-+-y-+4 i=0 
x—2 %—3 J 
Mee on as Ye | a Wok re — 
aS GA) t+y, vehge oa Vi ye i ny Yt+y,x 
2v—1 bs 1 it (94) 
a YS: be ee Lae VER pe 
2v (xX—1) p—3 2v (x—1) B ai 
% 2yv—] t+yv—1 
> f Xx Yirrit SY @WYie 
i=t+y+1 i—0 
FOR) 4.2, 
(2v—1)(x—2 eo 2 aN t--v—| 
Gi —Vi eee) EAN gee ae es as ba Yi,x—1 + >> (ft) Yix|— 
Qv (x—1) 7 Koren Yh | Fei dy i—0 
Ee wa ae ag ay Lea es ate 


2v—1 


ie ie (95) 
Dy (x—1) a (x) ae: R Dv ( on b (r--4, x) Yt+y,z i 


% 2v—1 t+v—1 
+ % «| ey tL k—te > o¥e| 


R==2 1=t +y+1 xX a EO 


for 2 x~<%, and at last for x =x, we have a formulae analogous to (94). 
We shall denote by X,a? and by Xo? the sum of coefficients squared 
of (94) and (95) respectively. We have 


Si eeu) 1.082) 8) 8 (96) 
2v (x—1) 16v? (x—1)? 
ge A a alo cl (97) 


2v (x—1) 16v? (x—1)? 


Dividing respective (Ux,i— Vz)? into 2 ae when x=2 0r x=% or 
into peel in other cases and taking the mean of these ratios we get the 
estimate of s*. We see however that the difference 


1 
= Se apaeee 98 
Ba? — By? — (98) 
is not a very important one. Therefore we may calculate the estimate 
of variance co? also in calculating simply the variance of U’s and dividing 
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—_— 


it over %,o2. We must remember however that by using this method 
we risk to superappreciate s? a little. It is clear that the above formu- 
lae (96) and (97) are valid for {> v. To get the corresponding formulae 
for t=v we have only to put everywhere in (96) and (97) « instead 
of x—1. We shall have 


(%) 2 __ Qv—1) 1), @ = 1)? + — 2 8 


* h Qy % a 16v? x? (99) 
(vy) (2v—1) (x—1) (x — 2)? —8 

DY Gs fe TE as a 100 
3 2v % an 16v? x? p10) 


Now we proceed to the calculation of the estimate of variance 
of Vitp— Vt. As we have seen above the problem consists mainly in 
calculating })° the sum of coefficients squared of the expression Vitp—Vt. 
The calculations being simple and analogous to those above we omit 
them, giving only the finite results. 


Case Il O0<t~<i+p~<». Using the formula (93) writing the 
analogous formula for Vt+p and subtracting, we have 


* % 
Vitp — Vi=—— Y YHpyx-— YD Vtx— 
*—1 yo %—1 9 
1 
———~— —- (Yt pty x— Yisp+yice Yin — Fit) (101) 
Av (x—1) } 
1 yep pid 


Ta Ss  WYin— Yi) 
YO) ot 


The corresponding sum %2 will be: 


a 


2 2p—1 
Si) 2. ee oe oe 
Anal sc ‘ 4y2 (x—1)? (102) 
Clearly the same formula is valid for vO t<i+tp <2v—1 
Case Ila. tS v<t-+p; p<. In analogous way: 
2v—1 1 *~!1 1 
Villy — Ve Sg 1 OV x tt 
TP Oy (x—1) +P, meat ik +p,x a Oy (x—1) t+p.x 
1 19 es 2v—1 1 7 
Te Grn rm sms 
SF) ae ee tx Gay Ue Leer (Yt+p—y,1 
(103) 
| l+y—1 ' 
Yttp—yut Ytt+y,1—Yttyn) — Pie ps (t,t+p) (Vi,i —Yix) 
2v (1) pti yt 
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and, as the reader will easily find: 


Sy eee oe es, wr Der tay (104) 
*—1 yv(r—1)* 4 (x—1)? 


Case IIb t<v<t+p; p>y. We have 


x—1 1 % 
py Vidgets > Ytx — 


il (Some 1 4=? 


Vitp — Vit = 


%~— 


1 1 
Bree ay 
y 


a1) (Yity.1—Yttyx) — 


(vot) (105) 


y t+y—1 
— aaa py (Yi,i— Yi, *) 
2901) ptt py tt 


and 
ay ee 2 2 (2y—p)—1 


- (106) 
%*—1 4y? (x—1)? 


As it was in the case when the number of compared varieties 
was odd this Case IIb is in fact identical with the Case I. Let 
‘us introduce a new variete d the smallest distance between the plots 
sown with the varieties the ¢-+ p-th and the /-th. 


In the Case I and Case Ila 
d=p (107) 
and in the Case IIb 
d=2v—p (108) 
We obtain therefore the following result: 


If d<v and if the two nearest plots sown with the compared 
varieties are not separated by the v-th variety (Cases I and IIb). 


® 


a | 
" cau (109) 


%*—1 Av? (x—1)? 


Ly)? = 


If howerer the two nearest plots with compared varieties are 
separated by the v-th variety, the formula (104) is valid. As we can 
easily see, in this Case we have a better accuracy of comparison, %)? 
being less. 
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Case fe. o-—t— v3; p=% We have 


4 Vv 4y—1 + - 1 < y +e y 
t+y ae t a9 vee Yt+ ysl tf TVs ity, % on 
4y (x—1) = 1 —1) 
(110) 
ye: 4yv—1_ i ri 
Av (x—1) pee %*—'1 9 eae dy CaaS 2v (%—1) i ate Yen) 
q 2 2v—1 
SS Wikre (111) 


%—1 4yv*(x—1)* 


PAen 1s ase IIa. It remains to consider the case when one of the 
compared varieties is the ,-th one. Assume first 1<it~<v | 
We haveg 


! i yt FE yn ov Qvx—1 # 
v. i= : Yo ae se  ee XY Ygu— 
ye 2v% (x—1) ,.9 2vx(%—-1) 9 
. x—2 er 
ee VO Gee ape ae —_—_ 
oi 0,1 + Yo,x+1) DiGe: Te aembE ar Ag? 3c Ser Yt+y,2— 
lie 1 t+y—1 2y—1 
ae re eA Yn) + 
eee h=tty-+i (112) 
1 4 Qy—1 t+y—1 
“ Wes ps | +s YR, x—-1-- > t Yn 
21) oe ety 
Hence 
2 2v— Ay 2y?—Dxt1 
aes “Dea rae (118) 


para Qvu(u—1) Sy? ray (ee 


We leave to the reader the care to prove that the same formula 
is valid when v< 17. 

Case 'Illb. As the latest stage we shall consider the case when 
the y-th variety is compared with the first one corresponding to the 
subscript t=o. Omitting the: easy calculations, quite analogous te 
the preceding ones we write: | 


Ayu--4y—2%+2 Qv x—-Qv Ht ¥-1 
Vas ORV ie aE is a te ee ee a x EEE Op kc. 
y 0 4y% (x—4) GM Nios y,%) a Qvn (x1) 8 ver, 
sy, oy hee 1 a * iy ? ue y at oy! 
Sa), ee i, Pa Sis Ik 
are 0,1 e tt) 2v% (x—1) , 5 Pees ares es i= 
1 %. Zy—1 y—l : 
|. +————__ ¥ Dl (a et oe ae | (114) 
2y% (2-1) x=2 \i=y-+1 i" 
a 
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Considering this formula, one can easily see that it is different 
from (112) and we should expect that the corresponding “)? will be also 
different from (113). Nevertheless the reader will easily prove that 
in this case the sum )? has also the same expression (113). 

We have now got all necessary formulae. To calculate the results 
of an experiment by method of E. Zateski we have: 

1) To compute the level fol every x-th plot 


1 xty 
iy = enn 115 
evens i | bye) 
fae y 


when the number of compared varieties is odd 2v-+14 and 


1 X+y—1 
Lx a ee + CYakiy + Ya+ty) + ¥ Yi | (116) 
av i=x*—y-+1 


when it is even 2v, where Yi denotes simply the yield of the i—th plot. 

2) We compute further the excesses of every variety over the 
level. Let us denote again the yields of the ¢-th variety by Yz.x, x being 
the order of the plot sown by this variety. The corresponding level 
will be denoted by Zt,x. The excesses of the ¢-th variety will be given by 


vi,x = Yt,x — Lt,x (117) 


> 


_ 38) We calculate the mean excess of the ?¢-th variety 


= 
awh 
Le 


1 
Vieja), Vix (118) 


ms 


where m is the number of plots sown by the ¢-th variety for which it 
was possible to calculate the level and the excess vi,x. The sum ¥ is 
extended over all these excesses. x 

4) To compare the two varieties the s-th and the ?-th we com- 
pute the difference Vs — Vt 

5) It remains to appreciate the accuracy of the comparison. We 
start with computing the estimate of the variance of Y’s. If the num- 
ber of compared varieties is odd n=2v-+1 we have to distinguish only 
two cases, when the number of plots for which can be calculated the 
level is »—1 or « Taking the sums of corresponding to these two 


1 
cases (Ut,x— Vz )?, multiplying them into appropriate Vy? which may be ta- 


ken from the Table I, adding and dividing over (2v—1) (x—1)-++2z-the 
number of all summands, we get the estimate of the variance so’. 

If the number of compared varieties is even n= 2y then calcula- 
tions are more complicated, as to different yields of the same variety 
there correspond different M2. Besides there are two different values 
of Ya2 corresponding to the ’-th variety. We get the estimate of s? cal- 


27 
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(Ut x—V1)? 
a. 
taken from the Table V. 

6) Having the estimate of s* we proceed to the calculation of 
the estimate of the standard deviation of Vs — Vz. For this we have 
only to identify the considered case with one of several cases in our 
theory, to which there correspond different %)2. The value of corres- 


culating the mean of all , where the inverse of Ya? may be 


ponding y Zr? may be taken from one of our tables. Multiplying it into 
the estimate of s we get the estimate of the standard deviation of Vs — Vt. 
All this will be illustrated on two following examples. 


Examptie ‘I, The table represents the data 


of experiments with five varie- 


ties of beets. The first column 
gives the numbers of consecutive 
plots x. The second column con- 
A tains the percentages of sugar 
B | Y in beets of five different varie- 
Cy) 18,26 18,122 | + 0,188 ties A,B, C, 1, Z, sown on re- 
D| 17,93 18,114 | —0,184 snective plots. The third column 
EB} 18,07 18,092 | —0,022 — gives the levels Lx) calculated 
Al 18,23 18,038 | + 0,192 ; for all plots except for the two 
| B| 17,97 18,032 | —0,062 | first and tae two last ones, 
C| 17,99 18,004 | —0,014 | Computing the levels we fol- 
D| 17,90 17,930 | —0,030 | low the formulae given above. 
E| 17,98 17,934 | —0,004 § We start calculating the sums 
| ' of five consecutive numbers Y. . 
A| 17,86 17,848 | + 0,012} Having calculated the first sum, 
|B} 17,99 17,774 | + 0,216 — to have the next we subtract 
C | 17,56 17,704 | —0,144 } from the first sum the number 
Dt 63 17,668 | —0,1388 § Y4,1—=18,27 and add the number 
Ps (Do 17,622 | —0,042 | Ya,—18,23 and so on. 
A 172568 17,492 | >t 0188.1 18,27 
Bai tc t6 17,440 | + 0,320 | 28,08 
C| 16,91 | 17,196 | —0,286 | peer 
; ; 17,93 
DiVePl i iar 17,180 | + 0,090 | 18,07 
E | 16,36 17,024 | — 0,664 | apse 
‘ 90,61 
A| 17,60 16,894 | + 0,706 | 18.27 — 
B!| 16,98 16,806 | + 0,174 | 18,23 + 
Cl NG-26 16,814 | —0,554 § 
'D| 16,83 | 16,638 | + 0,192 | 90,57 
|B | 16,40 a eee 18,08 — 
17.9% 
Ah AGS mat es 
90,46 
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In this way we get the sums +0,61; 90,57; 99,46 etc. which is 
necessary to divide over 5 to get the numbers of the column L(x). 

The next stage is to calculate the excesses vz.x over the level. 
We do it subtracting the numbers of the column L(x) from these of the 
column Y. Now we proceed to the calculation of mean excesses Vi. | 

For that purpose we add the excesses vt,x corresponding to every 


variety and divide the sums over the number of summands. 


V4 UR Vo Up UR 
0,138 + 0,184 — 
0,192 0,062 — 0,014 — 0,030 — 0,022 — 
0,012 0,216 + 0,144 — 0,138 — 0,004 — 
0,188 0,320 +. 0,286 — 0,090 + 0,042 — 
0,706 0,174 + 0,554. — 0,192 + 0,664 — 
“1,098 0,648 0,860 — 0,070 — 0,732— 
1,098 0,648 —0,860 
4=— = 02745; Ve = =0,162, Vo =—~ =— 0,172; 
—0,07 —0,732 
Vo =~ =— 0014, Ve= 2 =— 0,188 


The next stage consists in calculsting of deviations vt..— Vt and 
their squares. Remembering that we have to calculate the mean of 
(vt.x — Vt)? 

da 
and for varieties C and D we arrange* the calculations in two columns 


and that Ma? has different values for varieties A, B and F 


Varieties v— V | (vu— V) Varieties 
| 
A — 0,0825 | 0,006806 | C 0,310 0,096100 
0,2625 | 0,068906 0,158 0,024964 
0,0865  0,007482 0,028 0,000784. 
04315 | 0,186192 — 0,114 0,012996 
| — 0,882 0,145924 
B 10.994 0,050176 
0,054 0,002916 D 70 0,028900 
0,158 0,024964 jj — 0,016 0,000256 
0,012 0,000144 — 0,124 0,015376 
- Pee) Seat 0,104 0,010816 
, ie a 0,206 0,042436 
0,179 0,032041 et a pees | 
0,141 | 0,019881 L(V it 
— 0,481 - |. 0,231861 = 0,378552 
|Z w—V)? = 
| = 0,656790 


MSE RATS Se NETRA Bit AEE EDIE TEI BOLE TEM OI NTS RSENS EE EMT PET ETE LER ET ENTS BAS LA LLL REIN NM OEE LF NET LY 
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We find the required inverses of sums Mg? in the Table I 
for v=2 and for two different values of x. lor the three varieties 
A, B, E,it is »=5 and for the others »—6. We have 


1 1 
Y, 2 = 18667 Y g2 = 1,625 


Thus the estim:ite of the varience 62 of Y’s will be 


0,656790 > 1,6667 ++ 0,378552  1,5625 


Esto? = — = 0,070648 
29 


and 
Esiz = Vy Esto? = 0,277 


Now we may proceed to the estimation of differences Va — Va 
ete. Considering the values of V’s we see that probably the best va- 
riety is A, the worse—#, but the significance of the respective diffe- 
rences can be estimated by means of their standard deviations. To get 
them we have to multiply the #sfg into a proper V D2. We start by 
noticing that in our case we shall not need the Table Hi as the 
distance between two plots separated by the plots sown with varieties 
the v-th and the v-++1—-st# is at least equal to 3, which in our case 
iSe= ye 2: 

Thus when comparing the varieties A, B and # among them we 
shall find the respective V}2 in the Table II, comparing one of 
these varieties with one of the remaining ones in the Table IV, 
and at last comparing the varieties C and D — in the Table II 
again, but assuming »=6, d= 1. 

Denoting by Rag the estimate of the standard deviation of the 
difference Va — Vg we have 


Rag. = Rap = Of711 4 sigs = 0 


The equality Raz = Rap is of course due to the fact that the 
nearest distance between plots sown with A and B and with A and # 
is the same: d=1. The same distance for varieties F and B is d=2, thus 


Rep = 0,714 X Ests = 0,198. 


and Ree is a little larger, but very slighty. Comparing the difference 
Va — Ve with Raz, we see that it is equal about 2,3 X Raz. We 
conclude that the difference is rather significant. So is the difference 
Ve — Ve, but the significance of V4 — Vg is doubtful — it may be 
easily the result of random sampling. Proceeding further we have 
Rac = Rap = Rec = Rap= Rec = RED = 0,663 X Estg == Oiiot 


The significance of the difference Vig — Vo = 0,417 = 2,42 X Rac 
39 


The Theoretical Basis of Different Methods of Testing Cereals 31 


seems to be sure. Very probably also is significant the difference 
Ve — Ve = 0,334 = 1,82 * Rac. The significance of other differences 
is doubtful. Now Rcp = 0,635 X Esto = (,176 and we see that the 
significance of Vp — Vc = 0,186 is also doubtful. 


Drawing IF 


29 ty @5 26 
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The drawings (IV) and (V) represent respectively the comparative 
height of yields and the level of fertility calculated as shown above. 
We see that even in that example with very unequal soil the applica- 
tion of E. Zateski’s method gives very good results. 

It. is worth noticing that the two autors E. Zateski and 
Mitscherlich consider the excesses vi,x as independent variates gi- 
ven directly by the experiments and appreciating the standard errors 
of differences Vs — Vi are using the mean value of (vt,x — V1)’ instead 

(vt,x — Vi y 
of our ee Ss Sera In this way they underestimate the standard errorr 


of Vs—Vt as Sx? is always < 1. For instance, in the considered 
ol 
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example the mean value of (vt,« — Vz)’ is 0,047061 and is considerably 
less than Esto? = 0,076643. 

Example Il. When considering the case with an even num- 
ber of compared varieties we have different formulae for x, 2? and dink? 
Nevertheless we omit the detailed description of way of calculations as 
it is very much the same as in the previous example. Let us compare 
four varieties A, B, C, D of sugar-beets sown on 21 plots. The three 
columns Y, Z and V are giving the observed yields, the levels and the 


excesses. 


Grouping tke excesses of dif- 
ferent varieties we easily calcu- 
late the mean-excesses 


2|B| 1818 Ve = — 0,072; Vp = — 0,185 
31 ¢C 17,94 18,06 BEF O13 The calculating of the estimate 
4|D 18,08 18,02 -L 0,06 of.g2 needs some few explications. 
. We group the squares (vi,x — Vt)* 
5| A 17,99 17,98 + 0,01 corresponding to different values 
6B) 17,984. 17,95. 6 60,08 [.. oftsiees 
710 | 7:80 LEST 0.07 
8|D| 17,97 17,78 | +°0,19 (vaz2 — Va)’ = 0,0068 
: (vA5 — Va)? = 0,0218 
9|A 17,50 17,69 — 0,19 (ve, — Vz)? = 0,0116 
10| B 17,58 17,63 — 0,05 § (vps — Vp)? = 0,0770 
pues e 17,63 17,62 + 0,01 (vpi — Vp)? = 0,0156 
12) Di 17,66 17,64 + 0,02 (vps — Vp)? = 0,0812 
18) iA | -17,70 17,50/94 92 0.114] YX, = 0,2140 
14; B| 17,60 17,54 | + 0,06 | 
gL La Sls eata Wea by 17,44, 1° 0273 (vas — Va)? = 0,0138 - 
16|.D| 17,70 1¥28° errOAT (vas — Va)? = 0,0338 
; (ves — Vp)? = 0,0008 
17) A} 16,90 | 17,11 | — 0,21 | (vn. — Va)? = 0.0189 
[| sere | imme | 088] om — Vo = 0 
; X, = 0,0945 
21; A| 16,36 | 
Inet wrnenecsceecmncern haere (vce — Vc)? = 0,0000 
(vc1 — Vc)* = 0,0023 (vcs — Vc)? = 0,0067 
(vcs — Vc)? = 0,0262 (vca — Vpay* = 0,0392 


Z, = 0,0285 XL, = 0,0459 


We have now to multiply %,, and X, into “ 


and X, and Xj 
cl | Og * 


remembering only that for the three varieties A, B and D the 
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corresponding value of « is x = 5, while for the variety C it is » = 6. 
Thus the estimate of the variance o? is 


,__ 0,214%1,7491+0,0945X1,7734+0,0285%1,6326+0,0459X1 6496 
0 aa a an a i 


and hence 
Ests = V Esto? = 0,198 


Proceeding to the estimation of standard errors Rag of differences 
Va — Vg we notice again that we shall not need the Table VII as 
vy being equal to 2 the compared varieties cannot be separated by the 
v-th variety. For the estimation of Rag and Rap we find the corres- 


ponding Y¥)2 in the Table VI forx=5, »=2, d=1. We have 


Ade = RAD OTE ESilg = 0.141 


Further-estimating Rac—=Rsc=Rcp we find the proper V ¥)2 in the 
Table IX 


Rac= Rac = Rep = 0,642 X Est 6 = 0,127 


At last, estimating Rap the distance between plots sown with 


varieties B and D being equal to v, we find the V S22 in the Table VIII 
and have 


~~ 


Rap = 0,699 X Est g = 0,138 


Considering the obtained values of R’s we see that we have no 
means of distinguishing among the three varieties A, B and C and that 
it is only possible to hope that some of the differences 


Vo— Va = 0,258; Vo— Ve=0,262; Vpo— Vo=0,257 


are significant. All these differences being nearly equal are about 
1,8 — 2 times their respective standard errors. If the frequency distri- 
bution of differences Va— Vp in repeated samples were normal, (which 
is very probable) then the probability of having such or less probable 
results, calculated under the assumption that there is no significant 
difference among all four varieties, would be about 0,03 — 0,07. This may 
be considered as sufficient reason for assuming that the fourth variety 
D is probably a better one that the three others A, 8B, C. 

We finish here the description of E. Zateski’s method noticing 
that the given formulae for estimating the standard errors are not the 
same as E. Zateski is using himself, so the eventual criticism should 


be addressed to the author. 
oo 
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What has been done in the present paper consists in the pointing 
out that the application of E. Zateski’s method assumes that 
the level of the soil-fertility satisfies a certain equation in finite diffe- 
rences. About this equation it has been shown that every straight line 
is its solution and that there are no solutions with distant maxima and 
minima. In the last part of the paper formulae are deduced giving 
the estimates of standard errors. 

We shall return to E. Zateski’s method in the third paper 
of the present series in which we shall compare the advantages of three 
different methods of calculation of agricultural experiments. 

The Tables for the present paper an computed by Miss K. Jaro- 
szynska. 
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Bol: Hy sls 
n=2v-+1 
1 (2v-+1) («—1) 
Se Me Ss (xO) 
4 2,2500 1,8750 1,7500 1,6875 
5 2,0000 1,6667 1,5556 1,5000 
6 1,8750 1,5625 1,4583 1,4062 
7 1,8000 1,5000 1,4000 1,3500 
8 1,7500 1,4583 1,3611 1,3125 
9 1,7142 1,4205 1,3333 1,2857 
10 1,6875. 1,4062 1,3125 1,2656 
11 1,6667 1,3880. 1,2962 1,2500 
12 1,6500 1,3750 1,28383 | 1,2375 
13 1,6363 1,3636 12724 1,2272 
14 1,6250 1,3541 1,2638 1,2187 
15 1,6153 1,3461 1,2564 —«-1,2115 
164 3 1,6071 1,3392 | 1,2500 ——-1,2053 
17 1,6000 1,3333 1,2444 1,2000 
18 1,5937 1,3281 1,2395 1,1953 
19 1,5882 1,3235 1,2352 1,1911 
20 1,5833 1,3194 1,2314 1,1875 
21 1,5789 1,3157 1,2280 1,1842 
Ne ce 


Do 
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n=2v+1 


the v-th and by the v-+1—st varieties. 


Bb 


vate ‘ 2 2d 
ho Vs T (eel) oe ae 
Mee yess n=2v+1=7 

ae 1 2 1 2 | 3 
4 | 0,8219| 0,8273 | 0,8186 | 0,8220 | 0,8247 | 0,8182 | 0,8198 | 0,8215 
5 | 0,7106 0,7142 | 0,7089 | 0,7089 | 0,7125 | 0,7082 | 0,7084 | 0,7095 
6 0,6350 | 0,6375 | 0,6337 | 0,6350 | 0,6356 | 0,6332 | 0,6340 | 0,6347 
7 | 0,5792| 0,5811 | 0,5783 | 0,5793 | 0,5808 | 0,5779 | 0,5785 | 0,5791 
8 | 0,5360| 0,5375 | 0,5354 | 0,5360 | 0,5368 | 0,5350 | 0,5355 | 0,5360 
9 | 0.5012} 0,5025 | 0,5006 | 0,5011 | 0,5019 | 0,5003 | 0,5007 | 0,5011 
10 0,4724 | 0,4734 | 0,4719 | 0,4724 | 0,4730] 0,4717 | 0,4720 | 0,4723 
11 | 0,481! .0,4490 | 0,4477 | 0,4481 | 0,4486 | 0.4474 | 0,4477 | 0,4480 
12 | 0,4271| 0,4279 | 0,4267 | 0,4272 | 0,4276 | 0,4266 | 0,4268 0,4270 
13 | 0,4089) 0,4095 | 0,4085 | 0,4089 | 0,4093 | 0,4082 0,4086  0,4089 
14 | 0,3928/ 0,3933| 0,3924  0,3928/ 0,3931 | 0,3924 0,3925 | 0,3928 
15 | 0,3784| 0,3789| 0,3782 | 0,3784 | 0,3787 | 0,3780 | 0,3781  0,3783 

t 16 | 0,8657) 0,3681 10,354 | 0,3655 | 0,3658 0,3652 | 0.3453 | 0,8655 
17 | 0,3540| v,3544 | 0,3537 | 0,3540 | 0,3541 | 0 3537 | 0,3533 0.3589 
18 | 0.9434} 0,3439 40,3431 0,3434 | 0,3487 | 0,3481 | 0,3431 | 0,3482 
19 | 0,3386| 0,3841 | 0,3385 | 0,3336 | 0,3237 | 0,3333 0,8353  0,8384 
20 | 0,3246 0.8251 | 03245 0,3247 0,3229| 0 8245 | 0,3245 | 0,3246 
21 0,3159 | 0,3156 0,3166 | 0,816710,5161 | 


0,3164 


0,3161 | 0,3162 | 0,3162 5 
| i 


Case I d<»v: the two compared vurietes arc not separated by 


0,8232 | 
0,7114 
0,6350 
05797} 
0,5364 
0,5015 f 
0,4727 | 
0,4482 § 
0,4273 § 
0,4091 | 
0,3930 Ff 
0,3785 
0,3657 | 
0,3540 
0,3433 f 
0,3334 
0.32464 
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he vou eel Ls 


n=2v+1 


Case II The two compared varieties are separated by the v-th 
-and_the yv-+1—st varieties. 


/3 = V2Q@v+1)? @—1)— Qy+i—p 


(2¥-+1) (x—1). 


4 0,8023 0,8063 0,8055 
5 0,6978 0,7006 0,7000 
6 0,6260 0,6277 0,6273 
7 0,5723 0,5737 0,5737 
8 0,5306 0,5317 0,5315 
9 0,4969 _0,4976 0,4975 
10 0,4687 0,4695 0,4694 
11 0,4449 0,4455 0,4454 
12 0,4244 0,4249 0,4248 
13 0,4065 0,4069 0,4068 
14 0,3907 0,3944 0,3948 
15 0,3766 0,3769 0,3768 
16 0,3633 0,3642 | 0, 8641 
17 0,3525 0,3527 0,3526 
18 0,3420 0,3422 0,3421 
19 0,3324 0,3326 0,3325 
20 0,3236 0,3238 0,3237 | 


oT 


o8 


a8 


TAR Lt av 
Case IIl n==2yv-+1 


Ww 


yr 


0,3280 


ioe) 


One of the compared varieties is the v-th or the v+1 st one.. 
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Tevet) GOs 


0,4160 
0,992 


0,2844 


0,3707 


0,9586 
0,3475 
0,3375 
0,3283 


60,3198 


(Qy+1)%(@—1). 


0,757 
0,6667 


0,6025 


0,5589 


0,5156 
0.4842. 
0,4581 
0,4358 
0,4164 
0,3994 | 
0,3844 


0,3709 


0,3588 


0,347 
0,3377 
0,3284 


0,3199 
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4 | 2,8835) 1,9864 


5 | 2,4154) 1,7491 


6 | 2,1978) 1,6326 
7 |2,0709) 1,5590 


8 | 1,9884| 1,5119 


10 | 1,8935! 1,4551 
11 | 1,8563) 1,4345 
13, | 1,8318} 1,4206 


{ 


96) 1,4072! 


9 | 41,9842! 1,4792) 1,: 


26) 1,2621) 


83) 1,2570| 


ToAs 8.0 Epy. 


faa ae 


3 4 
1,7969] 1,7126 
1,5885) 1,5186 


| 


90 | 1,7288) 1,3578' 


‘| 2.3256 


1.2481) 1,1963 


83,2713 


2.6123 


2.1640 


2.0614 


ire 
i<e) 
CO 
i 
feeds 


41,9444 


41,9004 


4 9 2 OQ) 
1,828 | 


| 14,8054) 


1 6496 1,4928 
15718) 41,4306 
| 
1,5221) 1,8884 
| 

| 
1,4624' 1,3370 


1,4409, 1,3186 


179) 1,3609} 1,2500 


1,4879) 1,3585) : 
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LT AVB EL Se ev ir. 


nim 2 WW 
Case I d<vy. The compared varieties are not serarated by 


V8v2(4—1) + 2d—1 
2v (*—1) 


the v-th variety. /y)2 = 


4 0,8207 0,8183 0,8222 0,8175 0,8191 0,8216 


5 0,7100 0,7083 0,7108 | 0,7078 | 0,7090 0,7106 


6 0,6346 0,6333 0,6350 0,6330 0,6340 0,6350 


7 0,5787 0,5780 0,5794 | 0,5777 0,5785 0,5708 


8 | 0,5357 | 0,5350 | 0,5361 | 0,5348 | 0,5353 | 0,5360 
9 | 0,5009 | 0,5005 0,5012 | 0,50083 | 0,5007 | 0,5012 


0,4716 0,4719 0,4723 


0,4266 | 0,4268 | 0,4271 
0,4083 | 0,4086 | 0,4089 
0,3926 | 0,3929 | 0,3927 


0,3780 0,8777 0,3784 


16 0,3655 | 0,3631 0,3656 | 0,3652 | 0,8654 | 0,3656 


17 | 0,3539 0,3537 0,3540 | 0,3536 | 0,3537 | 0,3500 


18 | 0,8432 0,3431 0,3434 | 0,34830 | 0,8440 | 0,8484 
19 | 0,3386 0,3333 | 0,3887 0,3334 | 0,3335 | 0,3833 


20 0,3247 0,3245 0,3248 | 0,3246 | 0,8246 | 0,8248 


11 9,4480 0,4475 0,4481 0,4473 0,4477 | 0,4481 
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WA) BO Be VEE 
na 2 


Gane (la. 2 < py 


The two compared varieties are seperated by the v-th variety 


Vxje = V8" G4) 2p 1 


2y¥ (%*—1). 


4 0,7838 0,7937 0,7916 

5 0,6858 0,6925 0,6909 

6 0,6173 0,6220 0,6210 

7 0,5658 0,5693 0,5685 

8 0,5254 0,5282 0,5276 

9 0,4925 0,4948 0,4943 

10 0,4651 0,4670 0,4663 

11 0,4418 0,4435 0,4431 

12 0,4218 0,4231 0,4229 

13 0,4041 0,4054 0,4052 
| 14 0,3885 0,3897 0,3895 
15 0,3747 0,3757 0,3755 

16 0,3622 0,3631 0,3629 

17 0,3509 0,3517 0,3515 

18 0,3405 0,3418 0,3411 

19 0,3311 0,3318 0,3316 

20 0,3224 0,3230 0,3228 


41 


42 


Case Il-e 


TABLE VII. 


2 = 


Noe By 


ysv (x—1)—2v+1 


0,8036 
0,6987 
0,6270 
0,5729 


0,5810 


— 0,4971 


0,4688 
0,4450 
0,4245 
0,4066 
0,3907 
0,3766 
0,3640 
0,3525 


0,3421 


~ 0,3825 


0,3253 


2¥(%—1). 


0,8072 
0,7008 
0,6280 
0,5738 
0,5819 
- -0,4979 
0,4696 
0,4456 
0,4250 
0,4071 
0,3911 
0,3770 
0,3643 
0,3528 
0,3422 


0,3326 


0,8091 
0,7022 
0,6290 
0,5748 
 0,5825 
90,4982 
0,4700 
0,4460 
0,4253 


0,4072 


0,3914 


0,37 72 
-0,8645 
0,3529 
03424 
0,3328 


0,3239 


-— 
= 


ry 2 
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a or 


Case II. One of the compared varieties is the v-th ane. 


sere Vi ff2x(@e—1) @vx— 40 Eee hen 1S 


Te eo 8 


5 0,6420 | 0,6590 0,6685 0,6656 

6 0,5843 | 0,5968 0,6001 0,6014 

% | 05400 -0,5495 —0,5521 0,5533 

8 0,5043 0,5120 0,5142 0,5149 

: 9 0,4751 0,4814 » 0,4831 0,4839 
§ | 10 | \0,4502 0,455 ‘ 0,4571 0,457 
4 | | 11 7 — 0,4290 0,4336 0,4348 0,4854 
" . 1 ie 04105 | 04145 | 0,457 0,4164 
| “18 0,3942 0,8977 0,988 0,3991 
14 | 0,8797 0,3829 0,3837 0,384.2 

15 | 0,3667 0,3696 0,3708 0,3707 

% 16 0,3550 03576 | 0,582 0,3583 
i 17 0,3443 0,3466 0,3472 0,3475 
4 . | 18 0,3345 0,3866 0,3372 0,3878 
19 0,3255° 0,3275 |  0,3280 0,3285 

03190 |  0,3196 0,3198 
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Teoretyczne podstawy roznych metod obliczania doswiad- 
czen rolniczych. 


I Metoda E. Zateskiego. 


Niniejszy artykut jest pierwszym z trzech poSwieconych tej same} 
kwestji. Wyniki dajq sie strescié w sp.sd6b nastepujacy: 

1) Warunkiem koniecznym i wystarezajacym do tego, by metoda 
obliczania E. Zateskiego nie wprowadzala btedédw systematycz- 
nych, jest spefnienie przez poziom wydajnosci gleby réwnania rézni- 
cowego 

+y 
C2y el)” Yim a Ea) 


i=—y 


gdy liczba poré6wnywanych odmian jest n= 2v711 
hal 


i Ay Yx — Yx—y -{- Yuty +2 >y, Yui (2) 
i= —y+] 


Sayan — 29) 

2) Sréd rozwigzan ré6wnan (1) i (2) zawsze znajduje sie prosta. 

3) Rozwiazania (1) i (2) nie posiadaja odlegtych extremow. 

4) Jesli poziom wydajnosci gleby spelnia warunek (1), wowezas 
doktadnosé poréwnania odmian moze byé oszacowana na podstawie 
wzorow (70) do (118). 

Artykut nastepny bedzie poSwiecony metodzie obliczenia doswiad- 
ezen polowych, ktéra nazwiemy metoda krzywych parabolicznych. 

Artykul trzeci zawieraé bedzie poré6wnanie znanej metody ,Stu- 
dent’a“, metody E. Zateskiego i metody krzywych parabo- 
licznych. 
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